Above-Bandgap Magneto-optical Kerr Effect in Ferromagnetic GaMnAs 
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We have performed a systematic magneto-optical Kerr spectroscopy study of GaMnAs with vary- 
ing Mn densities as a function of temperature, magnetic field, and photon energy. Unlike previous 
studies, the magnetization easy axis was perpendicular to the sample surface, allowing us to take 
remanent polar Kerr spectra in the absence of an external magnetic field. The remanent Kerr angle 
strongly depended on the photon energy, exhibiting a large positive peak at ~ 1.7 eV. This peak in- 
creased in intensity and blue-shifted with Mn doping and further blue-shifted with annealing. Using 
a 30-band k ■ p model with antiferromagnetic s,p-d exchange interaction, we calculated the dielectric 
tensor of GaMnAs in the interband transition region, assuming that our samples are in the metallic 
regime and the impurity band has merged with the valence band. We successfully reproduced the 
observed spectra without ad hoc introduction of the optical transitions originated from impurity 
states in the band gap. These results lead us to conclude that above-bandgap magneto-optical Kerr 
rotation in ferromagnetic GaMnAs is predominantly determined by interband transitions between 
the conduction and valence bands. 

PACS numbers: 78.67.Ch,71.35. Ji,78.55.-m 



I. INTRODUCTION 

Dilute magnetic semiconductors ji i.e., II- VI or III-V 
semiconductors doped with transition metal ions, show 
a variety of magnetic phenomena arising from the in- 
terplay between localized and delocalized carriers. The 
strong s,p-d exchange interaction between band carriers 
(s,p) and Mn local moments (d) results in enormous en- 
hancement of carrier g-factors as well as the formation 
of magnetic polarons. The discovery of carrier-induced 
ferromagnetism in (III,Mn)V semiconductors^— with rel- 
atively high Curie temperatures (T c < 190 K)£ has stim- 
ulated further interest in these systems both from applied 
and fundamental viewpoints. 7 Although it has been ex- 
perimentally established that the ferromagnetic interac- 
tion between Mn moments is mediated by free holes^ 
there remain some basic questions as to the nature of the 
free holes (p-like or d-like) and the role of Mn impurity 
bands in transport and optical processes. 

Magneto-optical (MO) spectroscopy is useful for study- 
ing spin-dependent electronic states in magnetic systems 
via polarization-dependent reflection and absorption, i.e., 
the MO Kerr effect (MOKE) and magnetic circular 
dichroism (MCD). Previous MO studies of GaMnAs — 
however, have produced much controversy regarding the 
value and sign of the p-d exchange coupling constant, J p d, 
the position of the Fermi level, the insulating vs. metallic 
nature of the samples, and the contribution of the optical 
transitions originated from Mn impurity states. 

One of the commonly observed characteristics in MO 
spectra for ferromagnetic GaMnAs is that the sign of MO 
signal above the band gap (~ 1.4 to ~ 2 eV) corresponds 
to greater absorption for o~ polarization of light than for 



<j + polarization, which is opposite to (II,Mn)VI systems- 
as well as paramagnetic GaMnAs with low Mn doping. 
An opposite sign of MO signal could indicate an oppo- 
site sign of exchange coupling. Szczytko et al*^ proposed 
that the sign difference is a result of Ep being inside 
the valence band in ferromagnetic GaMnAs while p-d ex- 
change coupling is still antiferromagnetic (J p d > 0) as in 
(II,Mn)VI; due to exchange splitting, the Moss-Burstein 
shift 23 is expected to become spin-dependent, and thus, 
the lowest-energy a + transition should occur at a large 
momentum, while the lowest-energy a~ transition should 
still occur near the zone center. Komori et al. also ex- 
plained their MO data by adopting this models How- 
ever, Lang et al., 18 in analyzing their MO data within 
a parabolic band model assuming that Ep resides inside 
the valence band, had to conclude that p-d exchange cou- 
pling is ferromagnetic ( J pc i < 0); in order to explain their 
data quantitatively, they were also forced to introduce a 
dispersionless level of unclear origin inside the conduc- 
tion band that would need to have a very large oscilla- 
tor strength (three times larger than the usual valence- 
band-to-conduction-band transitions). In order to iden- 
tify the optical contribution from the impurity band, in- 
frared spectroscopy studies have been carried out as it 
probes electronic states near the Fermi surface. Burch 
et al . 16 ^ 17 interpreted the observed peak around 0.2 eV 
as the valence-band-to-impurity-band transition. On the 
contrary, Acbas et al. 22 measured infrared Kerr and Fara- 
day effects in GaMnAs and explained the spectra with a 
valence-band-hole theory. In addition, Ando et al£& and 
Berciu et al£^ both performed a MCD study on paramag- 
netic and low-T c ferromagnetic GaMnAs samples. Ando 
et al. claimed that the MO features below 1.4 eV are 
likely to be related to impurity states. They interpreted 
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their above-band-gap spectra in ferromagnetic samples 
as being composed of a broad positive impurity-related 
background and negative and positive peaks associated 
with band-edge singularities (-Bo and Eq + Aq critical 
points, respectively), concluding that E-p is in the Mn im- 
purity band and there is no Moss-Burstein shift. Berciu 
et at, based on their sample- dependent analysis and sim- 
plified one-dimensional model, proposed a unified inter- 
pretation of MCD spectra taking into account both im- 
purity band and valence band contributions. Different 
from Ando et al, they concluded that the broad posi- 
tive feature was due to valcnce-band-to-conduction-band 
transitions while the negative MCD near 1.4 eV, which 
occurred only when the sample was fully magnetized, was 
due to impurity-band-related transitions. 

Here, we report results of a systematic MOKE spec- 
troscopy study of GaMnAs samples with different Mn 
densities with three continuously varying experimental 
parameters: temperature (10 K - 75 K), magnetic field 
(-400 - +400 Gauss), and photon energy (1.4 eV - 
3.1 eV). Figure Q] displays representative Kerr rotation 
data measured for a Gao.976Mno.o24As sample (S-3) at 
10 K when the magnetic field was swept between positive 
400 Gauss and negative 400 Gauss. Strong variations of 
the Kerr angle are seen as a function of wavelength and 
magnetic field. Magnetic hysteresis is also evident by 
comparing the down-sweep [in (a)] and up-sweep [in (b)] 
data. 

We focus on the remanent Kerr angle, i.e., the Kerr 
angle at zero external magnetic field, which was found to 
strongly depend on the photon energy, exhibiting posi- 
tive peaks near 1.7 eV and 3 eV and a negative peak near 
2.5 eV. Without an external magnetic field, the MOKE 
spectra were caused by the interband transitions between 
the spontaneously spin-split bands. The ~ 1.7 eV peak 
increased in intensity and blue-shifted with the Mn den- 
sity and further blue-shifted with annealing. Using a 
30-band k • p model with antiferromagnetic sign of p-d 
exchange interaction, disorder broadening, and density 
dependent band gap, we calculated the dielectric tensor 
of GaMnAs in the interband transition region and suc- 
cessfully reproduced the observed spectra assuming that 
our samples are in the metallic regime and the Mn impu- 
rity band has merged with the valence band. There was 
no need to add any background dielectric constant or the 
optical transitions involving impurity states located in 
the band gap in order to explain our spectra. 



II. SAMPLES AND EXPERIMENTAL 
METHODS 

The three GaMnAs samples studied were grown by 
low-temperature molecular beam epitaxy (LT-MBE). 
They have a similar structure consisting of a 50-nm epi- 
layer of Gai-^Mn^As and a buffer layer of 1000-nm 
Ino.14Gao.86As on a GaAs (001) substrate. The sample 
structure is depicted in the inset of Fig. [2J Sample S-l 
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FIG. 1: Contour plot of the measured Kerr angle as a 
function of magnetic field and wavelength for an annealed 
Gao.976Mno.o24As (S-3) at 10 K. Magnetic field was swept 
from +400 Gauss to —400 Gauss in (a) and reversed in (b). 
(c) A hysteresis loop with magnetic field swept from +400 to 
—400 Gauss (solid line) and from —400 to +400 Gauss (dashed 
line), (d) and (e) show wavelength dependence of remanent 
Kerr angle (i.e., no external magnetic field) for down-sweep 
and up-sweep, respectively. 

had a nominal Mn content (x) of 0.01 and T c = 30 K. 
Sample S-2 had a nominal Mn content of 0.024 and T c 
= 45 K. The third sample, S-3, was from the same wafer 
as S-2, but was further annealed in air at 190°C for four 
hours. Its T c was increased to 70 K. Annealing is known 
to remove the interstitial Mn defects (Mn;) in LT-MBE 
GaMnAs and increase the effective Mn content x and the 
hole density p2^r— As a donor, Mni compensates accep- 
tors and decreases the hole carrier concentration in the 
sample. Mn; also prefers to be antiferromagnetically cou- 
pled to substitutional Mn su b in such a way that it reduces 
the local spins in the lattice. Thus, we expect that sample 
S-3 has a higher concentration of holes as well as effective 
Mn e ff than S-2. The sample parameters are summarized 
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FIG. 2: Temperature dependence of magnetization for sample 
S-2. A small magnetic field of 50 Oe was applied perpendic- 
ular to the sample surface. T c — 45 K is revealed. Inset: 
sample structure. 
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FIG. 3: X-ray diffraction patterns for three GaMnAs samples: 
(a) S-1, (b) S-2, and (c) S-3. Peaks due to the diffraction from 
the (004) planes are measured. InGaAs and GaAs peaks are 
observed around 65.2 and 66 degrees for all three samples. 
However, the GaMnAs peak is only observed in S-1 around 
66.5 degrees. 



sample 


S-1 


S-2 


S-3 


X 


0.01 


0.024 


0.024 


Tc(K) 


30 


45 


70 




un annealed 


unannealed 


annealed 



TABLE I: Gai-^Mn^As samples studied in this paper. S-3 is 
annealed and the other two are not. 



in Table UJ All the samples exhibited an "out-of-plane" 
easy axis, i.e., the magnetization direction was perpen- 
dicular to the sample surface. Temperature dependent 
magnetization of sample S-2 measured with a small per- 
pendicular field of 50 Oe is shown in Fig. [5J A T c of 45 K 
is clearly observed. 

The "out-of-plane" anisotropy indicates that all three 
samples are strained. The x-ray diffraction (XRD) pat- 
terns of the three samples are shown in Fig. [31 The 
spacings of the (004) planes of the InGaAs, GaAs, and 
GaMnAs layers are measured. In the case of sample S-1, 
the peaks of the three layers are observed at 65.2, 66, 
and 66.5 degrees, respectively, indicating that the lattice 
parameter of InGaAs is larger than that of the GaAs, 
whereas the lattice parameter of the GaMnAs is smaller 
than that of the GaAs. Therefore, the epilayer of sample 
S-1 is fully strained. The indium composition inferred 
from the lattice parameter is about 9.44% in the buffer 
layer. In sample S-2 and S-3, peaks for InGaAs and GaAs 
remain around 65.2 and 66 degrees, respectively. How- 
ever, the peak for GaMnAs is missing. It indicates that 
there is no well-defined lattice parameters for GaMnAs 
in these two samples. This happens when the epilayer is 
neither fully strained nor fully relaxed. Instead, the lat- 
tice parameter changes continuously across the growth 
direction. Thus, diffracted x-ray waves do not result in 
the well-defined Laue function with a sharp feature. 

MOKE measurements were performed in the polar ge- 
ometry in which the light beam was quasi-normally in- 
cident on the sample surface. Figure Ufa) is a schematic 
diagram of the experimental setup. White light from a 
100 W Xe lamp is first focused into a monochromator. 
Light with selected wavelength is then polarizaed with 
a Glan-Thompson polarizer and modulated with a me- 
chanical chopper. The linearly polarized light impinges 
on the magnetic sample at a nearly normal incident an- 
gle. The reflected beam passes through a photoelastic 
modulator (PEM), and an analyzer, and then its inten- 
sity is detected with a Si photodiode. The current signal 
is amplified and converted into a voltage and fed into 
lock-in amplifiers. The Kerr angle caused by the differ- 
ent reflection of is proportional to the signal at the 
second harmonic of the PEM modulation frequency. The 
two lock-in amplifiers are used to demodulate the signal. 
The first lock-in amplifier is referenced to the chopper 
frequency to provide a measurement of the average light 
intensity at each wavelength. The second lock-in am- 
plifier is referenced to the second harmonic of the PEM 
frequency to record the fast oscillating signal at 100 kHz. 
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FIG. 4: (Color online) (a) Schematic diagram of the exper- 
imental setup. Light with a selected wavelength is linearly 
polarized, and quasi-normal incident onto the sample. The 
reflected beam from the sample passes through a photoelastic 
modulator (PEM) and an analyzer, and it is detected by a 
photodiode. The signal is amplified and demodulated by two 
lock-in amplifiers. The sample is mounted inside a cryostat. 
Magnetic field is swept within ±400 Gauss, (b) Schematic di- 
agrams showing the directions of light propagation, magnetic 
field, and Kerr rotation. The incident beam propagates in the 
positive z-direction and is linearly polarized in the x-direction. 
The reflected beam travels in the negative ^-direction. It be- 
comes elliptically polarized and acquires a Kerr rotation 8k 
away from the ^-direction as shown in the right two panels 
in (b). 8k is defined to be positive (negative) if the rotation 
is clockwise (counterclockwise) when looking in the positive 
z-direction. 



The Kerr rotation angle (9k) is derived from the ratio of 
these two. 

Figure |4jb) schematically shows our configuration in 
more detail. Note that the magnetic field vector is oppo- 
site to the wave vector of the incident light. This is im- 
portant in determining how the refractive indices of circu- 
larly polarized normal waves in the material are related 
to the components of the dielectric tensor. The Kerr 
angle is defined as an arctangent of the negative ratio 
of the j/-component to the x-component of the reflected 
electric field, assuming that the incident wave propagat- 
ing in the positive z-direction is linearly polarized in the 
x-direction. To determine the sign of rotation, we first 
replace the sample with a silver mirror and induce a ro- 
tation of polarization by rotating the polarizer in front 
of the cryostat. Depending on the direction of rotation 
(clockwise or counterclockwise), there is a 7r-phase shift 
in the lock-in reading. We define it to be positive if the 
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FIG. 5: (color online). (a)Hysteresis loops (720 nm) detected 
for sample S-2 and S-3 at 10 K, 30 K and 43 K. Annealed 
sample S-3 exhibits larger saturation Kerr angle and coercive 
field at all temperatures. Temperature dependence of the (b) 
remanent Kerr angle and (c) coercivity for S-2 and S-3. 



rotation is clockwise. When the sample is mounted and a 
magnetic field is applied in the negative z-direction, the 
polarization rotation, i.e., the Kerr angle, is then mea- 
sured and recorded. 

Samples are kept in a helium- flow cryostat, allowing 
us to vary the temperature (T) from 10 K to above T c . 
An external magnetic field (H) is applied perpendicu- 
lar to the sample surface. The coercivity and saturation 
field are known to be small for all three samples from 
SQUID measurements. Thus, the magnetic field is swept 
within the range between ±400 Gauss. The Faraday ro- 
tation induced by the cryostat window was subtracted. 
Any polarization anisotropy caused by components in the 
setup was carefully calibrated and subtracted to get the 
accurate Kerr rotation. 



III. EXPERIMENTAL RESULTS 

Hysteresis loops measured with 720 nm (1.72 eV) for 
sample S-2 and S-3 are compared in Fig. 03a) at three 
temperatures, 10 K, 30 K, and 43 K. Generally, the an- 
nealed sample S-3 has a larger hysteresis loop than S- 
2 with both larger saturation Kerr angle and increased 
coercive field H c . At low temperature (10 K), the dif- 
ference between the two is relatively small; the rema- 
nent Kerr angles are 0.4° and 0.38°, respectively, and 
the coercive fields are ~ 170 Gauss and ~ 130 Gauss, re- 
spectively. With increasing temperature, the difference 
becomes more and more pronounced. The sharp vertical 
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FIG. 6: (color online). Temperature dependent hysteresis 
loops detected with different photon energies for sample S- 
2. Both the magnitude and sign of the remanent Kerr angle 
change with photon energy. For clarity, the hysteresis loops 
are vertically offset. 



switching of Kerr angle near the coercive field indicates 
the high quality of the sample with strong perpendicular 
magnetization for S-2. Annealed sample S-3 also presents 
strong anisotropy out of plane although the change over 
H c occurs at a slightly slower rate when compared to S- 
2. Figures [5jb) and[5jc) display the temperature depen- 
dence of remanent 9k (at 1.72 eV) and H c , respectively, 
for S-2 and S-3. It is clearly seen that T c increases as a 
result of annealing. 

The photon energy dependent hysteresis loops, 9k ver- 
sus H, for sample S-2 are presented in Fig. [5] For each 
photon energy, three hysteresis loops measured below T c 
(T = 10 K, 30 K, and 41 K) are displayed, and the size 
of the loop is seen to shrink with increasing temperature. 
At each temperature, the vertical size of the hysteresis 
loop (i.e., the remanent 9k) shows considerable variation 
with photon energy, while its horizontal width (i.e., the 
coercivity H c ) remains constant. As we can see, the hys- 
teresis loops are sharp and well-defined at most photon 
energies but disappear at 2.30 eV and 2.70 eV, although 
the temperature is below T c . In addition, the sign of 
remanent 9k depends on the photon energy. With the 
same sweeping direction, the hysteresis loop at 2.53 eV 
varies in the opposite way to that at 2.07 eV and 2.82 eV. 

Remanent Kerr spectra (9k at H = versus photon 
energy) for sample S-3 at various T (from 10 K to 75 K) 
are shown in Fig.[7][a). Above T c , there is no Kerr signal 
left. Below T c , each spectrum contains a pronounced 
positive peak at ~ 1.72 eV and a weaker positive peak 
at ~ 3 eV. The spectrum flips signs in the range of 2.4 - 
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FIG. 7: (color online), (a) Remanent (H — 0) Kerr spectra 
at different temperatures for S-3. (b) Remanent Kerr spectra 
at 20 K for S-l, S-2, and S-3. The main peak at 1.7 eV shifts 
to higher energy with increasing Mn content and annealing. 

2.7 eV. With increasing temperature, the main 1.72 eV 
peak remains in the same spectral position, while the 
magnitude decreases and finally disappears above 70 K. 
Spectra for S-l and S-2 were essentially the same as S- 
3, except that the main peaks occurred at slightly lower 
energies. Spectra for all three samples at 20 K are plotted 
together in Fig. [7{b) . We clearly see that the main peak 
shifts from 1.65 eV in S-l to 1.7 eV in S-2 and to 1.72 eV 
in S-3, due to increasing Mn content (S-l to S-2 and S-2 
to S-3) and annealing (S-2 to S-3). 

IV. THEORETICAL MODELING 

A. Band Structure Calculation 

Our theoretical model is based on the k • p 
method^Zr— and treats the antiferromagnetic exchange 
coupling between the hole spins in the host semiconduc- 
tor and the substitutional Mn 2+ spins within the mean- 
field approximation.— ~— Recently, this approach was 
used to explain observed MO effects in GaMnAs in the 
metallic regime, especially in the infrared spectral region 
near the Fermi level located inside the valence band. 22 
Here we model our experimental magneto-optical Kerr 
spectra in the interband region with a similar approach, 
but using a full band structure obtained with a 30-band 
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k • p method^ 9 - instead of the effective band structure cal- 
culation based on the 2-level model, which includes only 
Tec and rgv/r7V (here we follow the Td double group 
notation throughout the paper). We show that this ap- 
proach provides quantitatively and qualitatively reason- 
able results in reproducing the measured Kerr spectra 
and magnetization. 

To model MOKE spectra up to ~ 3 eV, the extension 
to the full band structure is critical, since the 8-band 
k • p method based on the 2-level model, which is of- 
ten employed to explain the electronic band structure of 
metallic GaMnAs, is valid only up to k ~ 0.1 A -1 mea- 
sured from the T edge, in both the conduction band and 
the valence band. Optical transitions with energy ~ 3 eV 
involve states with k- vectors higher than the above value, 
where the 2-level model breaks down. Furthermore, the 
dielectric tensor in the spectral range of interest obtains 
significant contributions from transitions with energies 
higher than 3 eV. Simulations based on the 8-band k • p 
model have to include an adjustable static dielectric con- 
stant added by hand. 

Furthermore, within the 2-level k ■ p model, one needs 
to specify the Luttinger parameters in order to take into 
account the interaction of the valence band with re- 
mote bands. The 30-band Hamiltonian does not con- 
tain the Luttinger parameters; instead, the coupling be- 
tween bands is described in terms of the matrix elements 
and band edges that are retrieved from the experiments. 
Once the band structure is determined with the 30-band 
method, one can calculate the Luttinger parameter s 29 ' 30 
for the subsequent use in simplified 8-band calculations. 
This is an important advantage of the 30-band k • p 
method as compared to the 8-band k • p method. There- 
fore, the modification of the valence band interaction 
with remote bands due to antiferromagnetic Mn-hole spin 
exchange coupling, strain, and Coulomb interaction ef- 
fects in diluted magnetic semiconductors can be simul- 
taneously taken into account within the 30-band k • p 
model. Also, the lack of inversion symmetry is taken 
into account through the matrix element between T^c 
and r 8V /r 7V . 

The extended k • p method is known to be valid up 
to 5 eV above and 6 eV below the top of the valence 
band, covering interband transitions of energies up to 
11 eV. It has been successfully applied to calculating 
effective masses and g factors for the group IV and 
III-V semiconductors! 29 ' 36 The general shortcomings of 
the method are also known; e.g., limited experimental 
data for higher conduction bands, the numerical difficulty 
in ensuring the continuity between U and K symmetry 
points. The accuracy of the GaAs band structure cal- 
culations with the 30-band method is discussed, e.g., in 
Ref.—. Here, we do not reproduce all the complicated ba- 
sis orbitals, energy level structures, band edge energies, 
and matrix elements for pure GaAs in a 30-band model, 
instead referring the reader to Refs. [2^ and[30l 



Our model Hamiltonian matrix is expressed as 

T -"0 ' n strain ' -"spin -ox 

#stai„ = E<*^ Q ^>5 ain W 

a,P 

fl£in-« = 5^(t|J(r-R J )Sx.s iC |i> 

~ Ji]N M nSMn(i\Ssc\j) 

where 

(i\D^\j) = - (PaP ^ ij +Vf. (2) 
m J 

Hq is a 30-band k ■ p matrix, H % 1 ■ is the strain Hamil- 
tonian matrix, e„| ain is the strain tensor, in which all 
off-diagonal elements are zero; a and /3 are coordinates, 
(i\D a P\j) consists of the linear combination of the de- 
formation potentials^ m is the bare electron mass, and 
-ffspin-ex describes the spin exchange interaction between 
the substitutional Mn magnetic impurity spins and the 
itinerant charge carrier spins in the host semiconduc- 
tor within the mean-field approximation. The indices 
i and j run over all 30-band basis orbitals. The an- 
tiferromagnetic spin exchange constant is not zero 
only when the bases correspond to Tec, Tsv, and T7v; 
Jij = 54 meV-nm 3 for Tg,v/Tjy, = —9 meV-nm 3 for 
T6C— and iV Mn = 4/a 3 3aMnAs , where a Ga MnAs is the 
lattice constant of GaMnAs. The strain effects are in- 
cluded only for the bases Tec, Tsv, and T?y. In this 
approximation, the T edges of Tgc, Tgy, and states 
are shifted by strain in the same way as those in the 
8-band model. 

We qualitatively took into account many-body 
Coulomb interactions through the phenomenological 
band gap narrowing (BGN). The hole-occupied exchange 
spin-split bands were assumed to be rigidly shifted 
by an amount ap 1 / 3 , where a = 2.6 x 10~ 8 cm-eV and p 
is in cm~ 3 . 39 The proportionality constant is compatible 
with that used in Ref. I22I . Several iterations were neces- 
sary to obtain self-consistent positions of the r§v edges 
and the Fermi level. Still, the resulting BGN shifts are 
only qualitative; rigorous calculations will be attempted 
elsewhere. The disorder effect was also phenomenologi- 
cally described as broadening of interband optical tran- 
sitions (~ 100 meV half-width at half-maximum) in the 
linear dielectric response function. 

The thermal fluctuations of Mn spin ordering at tem- 
perature T — 20 K, at which the experimental data in 
Fig. [7jb) was obtained, can be roughly estimated by com- 
paring the measured amplitudes of the remanent Kerr an- 
gle as a function of an external magnetic field at 10 K and 
20 K. The effect of thermal fluctuations was included in 
the calculation of temperature-dependent dielectric ten- 
sors and Kerr angle spectra. The temperature depen- 
dence of other physical parameters, for example, the band 
gap, the lattice constant, the strain tensor, and the hole 
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FIG. 8: (color online). Band structure of Gao.95Mno.05 As, 
calculated through the 30-band k ■ p method described in the 
text, taking into account the Mn-hole spin exchange interac- 
tion as well as a biaxial tensile strain of £^ am = +0.4%. The 
right panel shows a magnified band structure near k — 0. 

density were assumed to be negligible in this narrow tem- 
perature range. 

We assumed that our samples only have the epitax- 
ial biaxial tensile strain (see Section [XXJ> , which breaks 
the crystallographic cubic symmetry due to the lattice 
mismatch between the GaMnAs epilayer and the relaxed 
Gain As buffer layer. The same strain parameters as for 
GaAs 40 were used for the GaMnAs epilayer except for 
its lattice constant. Depending on the Mn (2 ~ 5%) and 
In (9.5%) contents, the strain tensor e|| ram (i = x, y, z) 
varies from 0.4 to 0.57% for the MnAs lattice constant 
5.98 Aj£ We assumed that our experimental post-growth 
procedure does not produce any additional significant 
strain. An example of the band structure calculated with 
the 30-band k • p method for 5% Mn is shown in Fig. [3] 



B. Hole Densities and Magnetic Anisotropy 

Before proceeding with the calculation of the 
frequency-dependent dielectric tensor and magneto- 
optical Kerr angle spectra, we roughly estimate the hole 
densities and Curie temperatures in our samples by cal- 
culating temperature-dependent magnetization based on 
the mean-field Zener model 3 - 3 .! 3 -^. In addition, we deter- 
mine the easy axis of magnetization with the estimated 
hole densities obtained for various Mn fractions. These 
will provide not only guidelines for the subsequent Kerr 
angle calculation but also justification for the results ob- 
tained by the extended kp method in the metallic regime 
in the interband range. To calculate the ferromagnetic 
critical temperature, temperature-dependent magnetiza- 
tion, and magnetic anisotropy, we use the 6-band k ■ p 
method since these quantities depend only on the very 





Mn (%) hole density, p (10 2U cm" 3 ) 
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2.35 


45 (S-2) 
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1.6 
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1.2 
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5 


1-7 r\ f c< n\ 

70 (n-3) 


4 


3.05 




5 


2.2 




2 


2.4 


30 (S-l) 


3 


1.39 




4 


0.97 



TABLE II: Hole densities calculated for various Mn fractions 
based on the 6-band k ■ p method and the mean-field Zener 
model for the measured Curie temperatures of our samples. 



top portion of the valence band structure. 

The temperature-dependent magnetization in Fig. [5] 
was calculated in the limit of strong degeneracy with zero 
external magnetic field based on the mean-field Zener 
model.™ Table \H\ lists the hole densities for various nom- 
inal Mn fractions, for which the calculated Curie tem- 
peratures are equal to the experimental ferromagnetic 
phase transition temperatures obtained from the SQUID 
measurement of the temperature-dependent magnetiza- 
tion in sample S-2 (Fig. [2]) and the measurement of the 
saturated remanent Kerr angle as a function of temper- 
ature (Fig. [5]). The calculation neglects the Fermi- liquid 
effect, which increases T c ^ as well as the spin- wave 
excitations^ which decrease T c for a given Mn fraction 
and a hole density. We simulate the annealing effect,- 25 
which increased T c to 70 K for sample S-3, by increasing 
the hole density for each Mn fraction. Note that generally 
not only the hole density increases, but also the lattice 
constant of GaMnAs is reduced upon annealing. 26 

The magneto-crystalline anisotropy, caused by the 
spin-orbit interaction, is attributed to the valence band 
holes since the total angular momentum of the local Mn 
moments is solely due to spins. The anisotropy gives cer- 
tain preferred directions for the spins to be aligned in 
crystals, so that the easy direction of magnetization is 
determined. We estimated the anisotropy field, which is 
proportional to the difference of hole free energies in [001] 
and [100] directions, 34 for different nominal Mn fractions, 
showing magnetic easy axes under tensile strain as a func- 
tion of hole density as shown in Fig. [TO] All hole densi- 
ties shown in Table HTl fall into the out-of-plane magnetic 
easy axis case ([001]) because they are all larger than 
the critical hole density for each Mn fraction, in which 
the easy axis is changed from [100] to [001]. The out-of- 
plane easy axis is further confirmed by the single domain 
like-behavior in the hysteresis loops of the measured re- 
manent Kerr angles for the applied out-of-plane external 
magnetic field shown in Fig. [3] 
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FIG. 9: (Color online) Calculated temperature-dependent 
normalized magnetization based on the mean-field Zener 
model and the 6-band kp method with different Mn fractions 
x = 0.03 (black solid line), 0.04 (red dashed line), 0.05 (blue 
dash-dotted line) for T c = 45 K (sample S-2) and 70 K (an- 
nealed sample, S-3), and with x — 0.02 (black solid line), 0.03 
(red dashed line), 0.04 (blue dash-dotted line) for T c = 30 K 
(sample S-l). The corresponding hole densities for each case 
are summarized in Table [Til The results are all for the ten- 
sile strain e^ rain = (0.57,0.52,0.46,0.4)% for Mn fraction 
x = (0.02,0.03,0.04,0.05). 



FIG. 10: (Color online) Anisotropy field, -ff an , calculated with 
the 6-band k • p model, as a function of hole density under 
tensile strain e^ rain = (0.57, 0.52, 0.46, 0.4)% for Mn fractions 
x — (0.02,0.03,0.04,0.05), respectively. The magnetic easy 
axes of our samples are out-of-plane ([001]) for the hole den- 
sities shown in Table [TTJ The curves are for Mn fractions 
x = 0.02 (dashed), 0.03 (dash-dotted), 0.04 (dotted), and 
0.05 (solid). 



C. Frequency-Dependent Dielectric Tensor and 
Kerr Angle Spectra 

With the estimated hole densities and magnetic easy 
axis directions, calculations of the dielectric tensor were 
performed based on the linear response theory (or the 
Kubo formula^—) shown in Eq. © for several different 
hole densities (shown in Table Mil from the range of the 
hole densities and Mn fractions calculated for our samples 
in the previous section (see Table HT) ) : 



=2fc2 



5a/3 



rngeo (27r) : 



£ 

a . b 



fa 



BZ 



(E ba y 



a P 
PabPba 



+ 



a P 
PbaPab 



(E ba - ih-y - frw) (E ba + ihj + hw) 



dfc a 



dkydk z 



(3) 



where a and (3 denote coordinates and a and b represent 
valence bands and conduction bands, respectively. Only 
interband transitions are included in calculations; f a is 
the Fermi-Dirac distribution and p" b are components of 
the momentum matrix element . The integration is ex- 
tended over the first Brillouin zone using full band struc- 
ture calculated with the 30-band k • p method described 
in Section |TV Al An example of the calculated dielectric 
tensor is shown in Fig. 111! 

Once frequency-dependent dielectric functions are ob- 
tained, calculating Kerr angle spectra is straightforward, 
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FIG. 11: (Color online) Complex dielectric functions e a 
e xy for 4% Mn, p = 1 x 10 20 cm" 3 and 5% Mn, p = 



-xy 
1Q 20 



x and 
1.5 x 



cm under tensile biaxial strain. 
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FIG. 12: (Color online) Calculated frequency-dependent 
magneto-optical Kerr angle under tensile strain by the 30- 
band k ■ p method in the interband transition range for vari- 
ous hole densities and three different Mn fractions iMn = 0.03 
(red dashed), 0.04 (black solid), 0.05 (blue dash-dotted) with 
GaMnAs epilayer thickness 70 nm. Hole densities (see Ta- 
ble IIIip were chosen with the range of ~ ±5 x 10 19 cm" 3 
around the hole density corresponding to T c = 45 K (sample 
S-2). For each Mn fraction, the higher peak corresponds to a 
lower hole density. 



taking into account all effects of propagation and reflec- 
tions in a thin GaMnAs layer and the buffer layer. For 
the dielectric function of the GalnAs buffer layer, the 
linear interpolation of the experimental data^ for GaAs 
and InAs was used. The calculated Kerr angle results are 
shown in Fig. [12] for sample S-2. Note that for a fixed 
Mn fraction, the Kerr angle peak magnitude decreases 
with increasing hole density. This happens because at 
these densities the Fermi level is already located below 
the lowest band edge of the spin-split Tgv- With in- 
creasing hole density, the Fermi level is shifted further 
downward, which reduces the number of states that con- 
tribute to interband optical transitions. On the other 
hand, if the Fermi level were initially in the middle of 
the spin-split Tgy, increasing the hole density would not 
necessarily decrease the Kerr angle amplitude. 

Kerr angle calculation results for 4% Mn give the best 
agreement with the experimental MOKE spectra of S-2, 
as shown in Fig. [12] The first positive peak that can 
be attributed to the interband transitions around the Eq 
critical point is red-shifted due to phenomenological band 
gap narrowing of 0.121 and 0.138 eV, respectively, for 
each case. Note that the above Mn fractions are 1.5-2 
times larger than the experimental nominal value, 2.4% 
for S-2, as we discuss in the following section. 

For sample S-l, the calculated result for Kerr angles is 
closest to the experimental MOKE spectra for 3% Mn, as 
shown in Fig. [13] and Table IIIII This Mn fraction is also 
larger than the experimental value of 1.5% for sample 
S-l. 

In Table HH the hole densities that are predicted for 
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TABLE III: Chosen Mn fractions and corresponding hole den- 
sities for T c = 45 K (sample S-2), T c = 30 K (S-l), and 
T c = 70 K (S-3) used in the dielectric function and Kerr an- 
gle calculations by the 30 band k • p method. The resulting 
Kerr angle spectra are shown in Fig. 1121 1131 and 1141 respec- 
tively. 
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FIG. 13: (Color online) Frequency-dependent magneto- 
optical Kerr angle under tensile strain calculated through by 
the 30-band k ■ p method in the interband transition range 
for various hole densities and three different Mn fractions 
x Mn = 0.02 (red dashed), 0.03 (black solid), 0.04 (blue dash- 
dotted) with GaMnAs epilayer thickness of 70 nm. Hole den- 
sities ( Table UTTj) were chosen around the estimated hole den- 
sities at T c = 30 K from Table [II] (sample S-l). For each 
Mn fraction, corresponding hole densities are smaller from 
the top. 
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FIG. 14: (Color online) Frequency-dependent magneto- 
optical Kerr angle calculated through the 30-band k ■ p 
method in the interband transition range for various hole 
densities and two different Mn fractions XMn = 0.04 (black 
solid), 0.05 (blue dash-dotted) with GaMnAs epilayer thick- 
ness of 70 nm. Hole densities (shown in Table HIT)) were chosen 
around the estimated hole densities (Table [UJ) for T c = 70 K 
(sample S-3) obtained by the 6-band k ■ p method and the 
mean-field Zener model. For each Mn fraction, corresponding 
hole densities are smaller from the top. 



T c = 70 K based on the 6-band k • p model are over- 
estimated since at large hole densities the hole free en- 
ergy calculated with the 6-band k • p model becomes non- 
negligibly smaller than the one calculated with the full- 
band model. This originates from the difference between 
the valence band structures calculated with these two 
models at large k. Therefore, it is expected that the full 
band structure calculation yields the hole densities lower 
than those listed in Table UT1 Note also that with 4% or 
5% Mn fraction, the Kerr angle calculation agrees better 
with the experimental Kerr spectra for T c = 45 K (S-2), 
as shown in Fig. [12j Therefore, for the annealed case of 
T c = 70 K the Kerr angle calculations are performed only 
for 4% and 5% Mn fractions. The results are shown in 
Fig. [TJJ and their corresponding hole densities are listed 
in Table |nll 



V. DISCUSSION 

Our calculated MOKE spectra of ferromagnetic GaM- 
nAs samples in the interband transition range reveal a 
general pattern of a large-amplitude positive peak around 
1.4-2.3 eV followed by negative and positive peaks with 
lower-amplitudes at higher photon energies, similar to 
the experimental data. At the same time, the quantita- 
tive spectral shapes and the positions of the peaks sen- 
sitively depend on the amount of substitutional Mn, the 
hole density, and the layer thickness. We determined the 
range of Mn fractions and corresponding hole densities 
for which the calculated Curie temperatures were equal 



FIG. 15: (Color online) Calculated MOKE spectra for ferro- 
magnetic GaMnAs with different layer thicknesses for 4% Mn 
fraction and hole density of 1 x 10 20 cm -3 . The dotted line 
shows the experimental spectrum for sample S-2. 



to the experimentally measured ferromagnetic transition 
temperatures. Then the parameters from this range pro- 
viding the best fit to the measured MOKE spectra were 
found. The resulting Mn fractions turned out to be 
larger than the experimental nominal doping values for 
all samples. If we assume that the experimental Mn con- 
centrations are quite accurate, this result could indicate 
that antiferromagnetic p-d exchange coupling strengths 
are 1.5-2 times stronger than the value of adopted 
in the simulations. This is because in the mean-field ap- 
proximation that we used throughout the calculations the 
exchange coupling Hamiltonian contains the product of 
2Vmh and Jy, as shown in Eq. (J]). 

Unlike magnetic circular dichroism, which has a clear 
physical interpretation as the difference of absorption co- 
efficients between the a~ and a + polarizations, it is dif- 
ficult to extract a single physical parameter that deter- 
mines the characteristics of the Kerr angle spectra in a 
thin layer limits since the contribution of the GaMnAs 
layer thickness, diagonal and off-diagonal components of 
dielectric functions, and the dielectric function of the 
buffer layer should be explicitly considered. However, 
the first positive peak can still be attributed to the in- 
terband transitions around the Eq critical point. The 
position of this peak and those of subsequent peaks are 
affected by the layer thickness. Figure [T3] shows a set of 
calculated MOKE spectra for a fixed Mn fraction (4%) 
and hole density (p — 1 x 10 20 cm~ 3 ) but different thick- 
nesses of the GaMnAs epilayer. The Kerr angle peaks 
red-shift as the thickness increases. The best match to 
the spectrum of sample S-2 is obtained when the layer 
thickness is about 70 nm, which is somewhat higher than 
our nominal value of 50 nm. 

For the annealed sample (S-3) the main peak in the 
measured MOKE spectrum shows a slightly increased 
amplitude as well as a slight blue-shift as compared to 
that before annealing (S-2). If we assume that the only 
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effect of annealing is an increased hole density for a given 
Mn fraction, our simulations predict an opposite trend: 
a decrease in the amplitude of the Kerr angle peak with 
increasing hole density. Within the mean field approx- 
imation, the observed effect of annealing indicates that 
the annealing also leads to an increase in the exchange 
interaction energy through an increase in the Mn fraction 
or /and J^-. 

VI. CONCLUSIONS 

We have used the magneto-optical Kerr effect to 
investigate three GaMnAs samples with out-of-plane 
anisotropy but different Mn densities. We measured the 
magneto-optical Kerr rotation as a function of three con- 
tinuously varying parameters: photon energy, magnetic 
field, and temperature. We observed remanent Kerr 
spectral peaks near 1.7 eV, which increased in intensity 
and blue-shifted with Mn-doping and further blue-shifted 
with annealing. We presented theoretical modeling of 
magnetic and optical properties of ferromagnetic GaM- 
nAs using the mean-field approximation and the 30-band 
k • p method, in which interband optical transitions are 



calculated over the whole first Brillouin zone and mod- 
ified valence band interactions with remote bands af- 
fected by the Mn-hole spin exchange interaction, the bi- 
axial strain, and the phenomenological Coulomb interac- 
tion and disorder effect are explicitly and simultaneously 
taken into account. Our calculation of MOKE spectra of 
thin-film ferromagnetic GaMnAs samples in the metallic 
regime is in good agreement with the experimental spec- 
tra. This supports the GaMnAs electronic band struc- 
ture in which the impurity band is smeared out into the 
valence band, although it cannot rule out the presence 
of impurity states in the band gap. To treat the MOKE 
spectra more rigorously in the interband range within 
the full band k • p method, the fc-dependent self-energy, 
which explicitly takes into account many-body Coulomb 
interaction and disorder effects, has to be calculated. 
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